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Intuitively one might expect independent noise to be a powerful tool for desynchronizing a population of 
synchronized neurons. We here show that, intriguingly, for oscillatory neural populations with adaptive 
synaptic weights governed by spike timing- dependent plasticity (STDP) the opposite is true. We found that 
the mean synaptic coupling in such systems increases dynamically in response to the increase of the noise 
intensity, and there is an optimal noise level, where the amount of synaptic coupling gets maximal in a 
resonance-like manner as found for the stochastic or coherence resonances, although the mechanism in our 
case is different. This constitutes a noise-induced self- organization of the synaptic connectivity, which 
effectively counteracts the desynchronizing impact of independent noise over a wide range of the noise 
intensity. Given the attempts to counteract neural synchrony underlying tinnitus with noisers and maskers, 
our results may be of clinical relevance. 

One of the important features of interacting neurons is the ability to coordinate their firing times and 
synchronize, which is essentially involved in normal brain function, ranging from sensory information 
processing to cognitive and motor function^'^. Abnormal neural synchrony, on the other hand, can 
significantly impair neural processes and is a hallmark of several neurological disorders^ Accordingly, syn- 
chronization control became a focus of the research nowadays^"^^. It is known that the extent of synchronization 
in oscillatory networks relates to the amount of coupling among oscillators^ An important mechanism that 
adjusts the synaptic coupling to the underlying neural activity is spike timing- dependent plasticity (STDP), where 
the synaptic weights are either potentiated or depressed depending on the order of the spike times of pre- and 
post-synaptic neurons^^"^^. This mechanism plays an important role in temporal coding of information by 
spikes^^'^^. 

Synchronized firing of neurons is characterized by narrowly distributed spike times, which, in turn, naturally 
leads to a potentiation of synaptic weights in neural ensembles with STDP and a stabilization of a strongly coupled 
and synchronized regime^'^^. On the other hand, if synchronization is perturbed, the synaptic weights are 
depressed until the oscillators get weakly coupled and thus desynchronized^'^^ This phenomenon motivated 
the development of stimulation methods for counteracting pathological neural synchronization. As shown in 
theoretical studies, animal experiments and in a clinical proof of concept study^'^^"^^, desynchronizing stimulation 
can lead to sustained desynchronization, pronounced normalization of effective connectivity and substantial 
decrease of symptoms. 

Desynchronization can effectively be achieved if the interacting oscillators with fixed coupling are stimulated 
by independent noise^^"^l It is thus to be expected that in neural ensembles with STDP independent noise will 
cause depression of synaptic weights. Indeed, in order to observe a synaptic potentiation due to STDP, a certain 
degree of correlation in the input is necessary^^. We here show that independent noise delivered to the neural 
population with STDP can surprisingly enhance the amount of synaptic coupling among the neurons, i.e., the 
neurons counteract the desynchronizing effect of noise by reorganizing their synaptic connectivity and increasing 
its strength. As a result, the extent of synchronization in a neural ensemble with STDP can be significantly larger 
as compared to without STDP if both receive an independent noisy input of the same intensity. Put otherwise, 
STDP leads to self- organized noise resistance. Since noise is inevitable in real systems and, hence, has an impact 
on neural information processing requiring precise spike timing^'^'^^'^^, our findings show how STDP may be 
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involved in a homeostatic mechanism preserving the existing level of 
collective dynamics in the brain by neutralizing the impact of ran- 
dom perturbations. 

Results 

We first demonstrate the discussed phenomenon for the ensemble of 
spiking Hodgkin-Huxley (HH) neurons (5) -(7) with STDP and ran- 
dom synaptic input (see Methods). If the input is absent {I = 0), the 
neurons get strongly coupled, and a uni-directional hierarchical 
coupling topology is established due to STDP for strong enough 
initial coupling [Fig. lb, Supplementary Fig. SI]. If an independent 
random input is applied to the HH ensemble, the amount of coupling 
in the perturbed ensemble is significantly enhanced, and, intrigu- 
ingly, the random perturbations have a constructive effect on the 
dynamics of the synaptic weights [Fig. la]. In fact, the noise pro- 
motes the development of bidirectional synaptic connections among 
neurons [Fig. Ic], the mean coupling K increases as the input 
strength grows, and there exists an optimal noise intensity, where 
the neurons get maximally coupled [Fig. la]. For the considered 
parameters the optimal input intensity I = Jopt ~ 0.14, where 
K~0.374, which corresponds to about 50% of relative increase of 
the mean coupling with respect to the input-free case. Furthermore, 
if the noise is not particularly strong, the amount of synchronization 
in the neural ensemble with STDP is relatively well preserved 
[Fig. la]. For example, for I = Jopt» the time-averaged synchroniza- 
tion order parameter (i^sTDp) ~ 0.S2 for the ensemble with STDP, 
whereas (i^pix) ~ 0.21 for the ensemble without STDP and for fixed 
coupling matrix Kpix (see Methods). This corresponds to about 11% 
and 77% of synchronization suppression relative to the input-free 
case, respectively. This mechanism constitutes a self- organized res- 
istance to independent noise, where STDP plays a stabilizing role for 
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the synchronized neural firing against the destructive influence of 
external random perturbations. 

The mechanism behind the reported phenomenon is based on the 
noise-induced dynamics of the spike time differences Aty used for the 
update of the synaptic weights [Supplementary Section SIB]. If 
the distribution density Pij{At) of Atfj is known, the average update 
rate of the synaptic weights can be calculated by the integral 



Akj = jw{0p,^{0d^, 



(1) 
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where W describes the update rule in one spike event (7). One there- 
fore has to describe the form of pij in order to predict the dynamics 
of the coupling weights kif. They will be potentiated if Akij>0 
and depressed if Akij<0. For instance, for a uniform distribu- 
tion = 1/28 of Atij e [ — 8,8], Akij = [AiTi{l— exp{ — 8/Ti)) — 
A2T2(1— exp( — e/T2))]/2£ and changes its sign from positive to 
negative as 8 increases over 5 ms for the considered parameters. 
Therefore, if Aty are broadly distributed, e.g., for desynchronized 
neurons at a weak initial coupling, the synaptic weights will be 
depressed and a weakly coupled regime will be realized. On the other 
hand, a synchronized firing for a strong initial coupling will lead to 
narrowly distributed Atfj and further stabilization of a strongly 
coupled regime [Fig. lb], see also Supplementary Section SI A. 
Notice that the considered STDP will lead to a potentiation of syn- 
aptic weights if the neurons fire at a high rate^^'^^. In this case, the 
relative spike times will be narrowly distributed even if the neurons 
are not synchronized. In our case the mean firing rate does not 
significantly change when the parameters of interest vary 
[Supplementary Fig. S3]. 

We consider the behavior of Py in more detail for the phase oscil- 
lators (8) with STDP, where we observe the same constructive effect 
of the independent noise: As the noise intensity increases, so does the 
mean coupling strength K until it reaches its maximal value at the 
optimal noise intensity /n = ^opt [Fig- 2a]. As for the HH ensemble, 
the noise induces the emergence of bidirectional couplings among 
oscillators [Fig. 2b], and STDP counteracts the suppression of syn- 
chronization by the noise [Fig. 2a]. Although the latter effect for the 
phase ensemble is somewhat less pronounced (regarding (i^sTDp)/ 
(i^Fix)) than for the HH neurons, we still observe the self- organized 
resistance to noise for the phase oscillators mediated by STDP, where 

(^STDP) > (^Fix)- 

We found that the emergence of positive couplings below the 
diagonal in the coupling matrix [Fig. 2b] of the phase ensemble 
can be explained by evaluating the impact of the noise on the phase 
differences for fixed coupling matrix Kpix, and we consider such a 
system below. 

In the strongly coupled regime without noise the oscillators are 
phase locked and the phase differences are constant. The following 
recurrent formula gives the constant phase shifts of the locked phases 
ij/i = wt + (pi (see Supplementary Section S2A): 



Figure 1 | Constructive effect of independent random input on the 
synaptic weights of the HH neural network with STDP. (a), Mean 
(ensemble-averaged) synaptic weight K (scale on the left vertical axis) and 
the time-averaged order parameters (-Rstdp) and (-Rpix) (scale on the right 
vertical axis) for the ensemble with and without STDP, respectively, versus 
input intensity I. For the network without STDP, the coupling matrix is 
fixed Kpix (see Methods), (b), (c), Coupling matrices established in the HH 
ensemble due to STDP for the input intensity 7 = 0 in (b) and I = 0.14 in 
(c) , where the synaptic weights kij are encoded in color. The elements of the 
initial coupling matrix K(0) = { A:y(0) } are Gaussian distributed around the 
mean value K(0) =0.5 with standard deviation 0.02. For illustration, the 
neurons are sorted with respect to increasing natural spiking frequency 
such that fi < fj for i < j. 



( Nwi 
+ arcsin - — , 



(2) 



where D,+ 1 = Sf^ ^ + Cf+ ^ , S, + 1 = Yl^^ ,+ 1 (^;) ' C/+ 1 = Ejl /+ 1 
cos (^(pj^ , and Wi = coi — con. Letting (p^^ = 0, without loss of general- 
ity, we find the constant phase differences Axj/ij = cpj — (pf. For weak to 
moderate noise the phase differences fluctuate around the above 
mean values and get broadly distributed, see Fig. 3a, Supplemen- 
tary Section S2B, and Refs. 25, 26. Since the phase differences obey 
a circular statistics we approximate them in the noisy case by the von 
Mises distribution^^ 
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Figure 2 | Constructive effect of independent noise on the coupling of the 
phase ensemble with STDP. (a), Mean coupling weight K and ratio 
(-RsTDp)/(-RFix) (scale on the right vertical axis) between the time -averaged 
order parameters for the ensemble with and without STDP, respectively, 
versus noise strength fi. In the latter case the coupling matrix is fixed Kpix 
(see Methods), (b). Coupling matrix Kstdp established in the phase 
ensemble due to STDP for the noise intensity fi = 0.1. The dashed black 
curve in (a) approximates the mean coupling weight K, and the white curve 
in (b) bounds the region of potentiated coupling weights, both as predicted 
by the theory (see text for details). The oscillators are sorted with respect to 
increasing natural frequencies Wj which are uniformly and randomly 
distributed in the interval [0.9, 1.1]. 



M{k,(p) = 



2kIo{k) 



exp(K:cos(^ — ^o))' 



(3) 



where Io(k) is the modified Bessel function of first kind and zero 
order, k is the concentration parameter, and (po is the mean value. 
The parameter K:(Ai/^y) can be determined and demonstrates a 
power-law dependence on the noise strength /n [Fig. 3b] 



COiCOj 

(hi + coj 



(4) 



where the exponent y ~ — 2 seems to be independent of the natural 
frequencies cOi and coj [Fig. 3b and Supplementary Section S2B]. We 
therefore represent the dependence of parameter K:(Ai/^y) on the 
natural frequencies as a function F(Wi, coj) in equation (4). Detailed 
numerical calculations reveal that the contour lines F(Wi, Wj) = F const 
of the function F align along the hyperbolas PwiWj — F const {^i + ^j) 
= 0 for some constant j], where P ~ —9.52 for the considered fixed 
other parameters of the system, see Supplementary Section S2C. 

With the use of equation (4) one can predict the noise-induced 
emergence of positive coupling from slow oscillators to fast ones, 
(for indices) < / in Fig. 2b). For this, the STDP function (7) can be 
integrated with respect to the distribution density p( Ai/^/y) of the form 
(3) according to equation (1) for the mean phase differences (2) and 
parameter k from equation (4) with the obtained parameters y and 
The results of such an approximation are illustrated in Figs. 3c, d, 
where the range of potentiated coupling weights increases as 
predicted by Akij>0 [Fig. 3c] which is in good accordance with 
the numerical simulations for weak to moderate noise [Fig. 3d]. 
Moreover, the above equations can effectively be used to approx- 
imate the mean coupling [Fig. 2a] and the shape of the coupling 
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Figure 3 | Impact of the noise-induced dynamics of the phase differences 
on coupling weights of the phase ensemble, (a). Noise-induced 
distribution densities pij of the phase differences Ai/^y(/i) between the post- 
synaptic oscillator i = 100 and pre-synaptic oscillators j indicated in the 
plot. The hatched histograms illustrate the results of numerical simulations 
for noise intensity ji = 0.1, the enveloping dashed black curves depict the 
von Mises distribution (3) for fitted k:, and the vertical lines are delta-peak 
distributions for the noise-free case according to equation (2). (b). Log-log 
plot of parameter k:( Ai/^y) of the fitted von Mises distribution (3) versus fi 
for oscillator pairs ( z, j) indicated in the plot. The dashed black lines have 
the slope —2. (c). The average update rate A/ciooj of the coupling weights. 
The black dashed curves depict the theoretical approximation calculated 
according to equations (l)-(4) with y = — 2 and P = —9.52, while the 
symbols represent the results of numerical simulations. Vertical dashed 
lines indicate the index j where the predicted A/ciooj crosses zero. In plots 
(a)-(c) coupling matrix is fixed to Kpix (see Methods), (d). Time-averaged 
coupling weights (/cioo,;) established in the phase ensemble with STDP for 
noise intensities as in plot (c). Natural frequencies cOi = 0.9 + 0.2(z — 1)/(N 
- 1), i= 1,2, ...,iV, N= 200. 

matrix with irregular border between potentiated and depressed 
coupling weights [Fig. 2b] for randomly distributed natural 
frequencies. 

Strong noise destroys the synchronization in the HH and phase 
ensembles, which, in turn, leads to a depression of the coupling 
among oscillators by STDP, and the mean coupling weight K nearly 
vanishes [Figs, la and 2a]. We found that the onset of decoupling 
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originates from the fast oscillators (large index j in Figs. Ic and 2b) 
and propagates to the interior of the coupling matrix as time evolves, 
see Supplementary Videos 1 and 2 for the animated time courses of 
the coupling matrices of the HH and phase ensembles, respectively. 
The indication of such a decoupling mechanism can be obtained 
from the population without STDP and fixed coupling matrix Kpix 
(see Methods). Already for /n = ^opt ^ 0.1625 the coupling weights 
from the fastest oscillators are predicted to be depressed by STDP as, 
for example, for / = 100 and; = 200, where AA;ioo,2oo <0 [Fig- 3c]. 
This indicates that the suppression of the coupling for large noise will 
start from the fastest oscillators. 

For the same parameters of the STDP function (7) the self-orga- 
nized resistance against noise is significantly enhanced for another, 
multiplicative (state-dependent) update rule, where the coupling 
weights are changed relative to their current values ky k^ + (/cm 
- kij)\SW{Ail/ij)l where /cm = k^^^ or /cm = k^^ depending on 
whether the coupling update dW{A\l/ij) is positive (potentiation) or 
negative (depression), respectively^^ Even for strong independent 
noise the oscillators with such a coupling update remain strongly 
coupled and synchronized [Fig. 4a]. The dynamics of the coupling 
matrix is, however, different for this update rule [Fig. 4b]. As the 
noise strength increases, the coupling weights k^ converge from a 
strictly bimodal distribution for the noise-free case via a nearly uni- 
form distribution to a well-pronounced uni-modal distribution of 
approximately globally coupled oscillators for strong noise. The lat- 
ter distribution is in accordance with the results of Ref 28. 
Interestingly, the oscillators are coupled with about the same mean 
coupling in the noise-free case (fi = 0) and for strong noise (fi = 0.5) 
[Fig. 4a]. In the latter case, however, the population with STDP is 
much more resistant to noise than the population without STDP, 
while both have the same mean coupling. This also indicates that the 
global coupling with narrowly distributed coupling weights [Fig. 4b] 
is more favorable if robustness of the synchronized dynamics against 
noise is required. These findings also contribute to the problem of the 
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Figure 4 | Constructive effect of independent noise on the coupling of the 
phase ensemble with multiplicative (state-dependent) coupling update. 

(a), Mean synaptic weight K (scale on the left vertical axis) and time- 
averaged order parameters (-Rstdp) and (-Rpix) (scale on the right vertical 
axis) for the ensemble with STDP and without STDP (for fixed coupling 
matrix Kpi^, see Methods), respectively, versus noise strength fi. (b), 
Normalized histograms of the distribution density of the coupling weights 
kij for the noise intensities = 0 (no input, scale on the left vertical axis), ji 
= 0.125 and ji = 0.5 (scale on the right vertical axis) as indicated in the 
plot. 



optimal coupling topology for easily synchronizable networks and 
their resistance to noise^^. 



Discussion 

We demonstrated the constructive impact of independent noise on 
ensembles of oscillatory HH neurons and phase oscillators with 
asymmetric STDP and for two different protocols of the coupling 
update based on the difference of either the spike timing or the 
phases. In the latter case the coupling update is separated from 
possible changes of the frequencies, which means that the reported 
phenomenon is in fact based on the phase dynamics. This mech- 
anism is different as compared to stochastic or coherence reso- 
nances^""^^, where the interplay between noise-induced and 
internal or external times scales is important. Indeed, the noisy input 
does not significantly change the frequency of the synchronized 
neurons and only influences the distribution of the phase differences 
by broadening it [Fig. 3] . This results in a potentiation of the synaptic 
weights that were depressed in the noise-free case. 

For the phase model we provided a theoretical approximation of 
the structure of the coupling matrix, which is in good agreement with 
the results of the numerical simulations. For the HH ensemble we 
also verified the reported phenomenon for a random synaptic inhib- 
itory input and somatic Gaussian noise as well as for different spike 
pairing for STDP [Supplementary Section SIC]. Furthermore, we 
verified the phenomenon for different parameter sets for the plas- 
ticity function and found that the noise-induced increase of synaptic 
weights gets more pronounced if the STDP function is shaped for a 
stronger potentiation rather than for depression. In this study we 
however considered a non -trivial case, where strongly coupled 
synchronized and weakly coupled desynchronized states stably coex- 
ist if the communication delay among neurons is small, and syn- 
chronization is associated with large coupling^^'^^. In general, the 
noise-induced increase of the mean coupling can be expected for 
an asymmetric STDP rule if the neurons can be in-phase locked with 
narrowly distributed spike time differences, and a uni-directional 
coupling topology like in Fig. lb is established in the noise-free case. 

The resonance-like behavior of the mean coupling versus noise 
intensity has been observed for both multiplicative and additive 
coupling update rules, where the coupling update either depends 
or not, respectively, on the current value of the synaptic weight 
and leads to "soft" or "hard" bounds, accordingly^^'^^. The distri- 
bution of synaptic weights however is always bimodal for hard 
bounds [Figs. 1 and 2], whereas it undergoes a gradual transforma- 
tion from a bimodal to a uniform distribution and then to a relatively 
narrow unimodal distribution as noise strength increases for soft 
bounds [Fig. 4], see also Refs. 23, 28. These findings point to the 
robustness and generality of the reported phenomenon. 

Our results suggest a possible homeostatic mechanism of how the 
brain may counteract external perturbations and noise in order to 
preserve the existing level of neural synchrony and bridge the trans- 
ition from information coding by precise spike times to variable and 
imperfect spike timing. Noise leads to the development of bidirec- 
tional coupling among neurons observed in experiments^^. Further- 
more, our results show that independent noise can by no means be 
considered as an effective method for desynchronization of oscillat- 
ory neural networks with STDP. Also from the clinical standpoint 
our results maybe important. In fact, they may contribute to a deeper 
understanding of why maskers and noisers show limited efficacy in 
counteracting tinnitus^^, the latter being associated with abnormal 
neural synchrony^'^°. 



Methods 

Hodgkin-huxley (HH) ensemble. The considered system of N synaptically 
interacting spiking HH neurons^^'^*^ reads 
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CVi = li -gmm'MVi - y^a) -gKn^AVr - Vk) 

_ (Vr-Vf)v",^_^^,^^. input 



-gi(Vi-Vi) + - 



N 



= am{Vi){l-mi)-l]^{Vi)mi, 

-a,{V,){l-h)-[^,{V,)h,, 

= 0in{Vi){l-ni)-l3^{Vi)ni, 

= 5(l-s,-)/(l + exp(-(y, + 3)/8))-s,, 



(3) 



where the variable Vi,i= 1, . . N, models the membrane potential of neuron /, and 
ajV) = (0.1 y+ 4)/[l - exp(-0.1 V - 4)], JV) = 4 exp((-y - 65)/18), a;,(y) 
= 0.07exp((-F- 65)/20),jS;,(V) = l/[l + exp(-0.1 V - 3.5)], a„(y) = (0.01 V + 
0.55)/[l - exp(-0.1 y- 5.5)],andjS„(V) = 0.125 exp((-y- 65)/80). Parameters C 
= 1 /xFW, V^a = 50 mV, Vk = -77 mV, Vi = -54A mV,^jva = 120 mSW,^^ 
= 36 mS/cm^^/ = 0.3 mS/cm^. 

The neurons are excitatorily coupled (with reversal potential Vr = 20 mV) via the 
post- synaptic potentials (PSP) with synaptic weights ky modeling the strength of the 
coupling from the pre-synaptic neuron j to the post-synaptic neuron i. We consider 
the case without self-coupling ku = 0. The constant currents are randomly and 
uniformly distributed in the interval [10.5; 11.5] fiA/cm^. For such parameters the 
spiking frequencies/ (number of spikes per second) of the uncoupled neurons {ky = 
0) without an external input (f^^^^ = o) get distributed around the mean value 
70.7 Hz with standard deviation 0.62 Hz. For numerical integration of model (5) with 
synaptic input we used a Runge-Kutta method of order 5(4) with adaptive step size^^. 



Independent random input. The excitatory synaptic input current if^^^ reads 



tr\t)=i{v,-v,{t))j2<t-c,,). 



(6) 



The a-train in equation (6) models the PSPs received by neuron i at times T,jt with 
inter-spike intervals At/^^ = Zi^k+i ~ ^i,k — 0 independently drawn from a Gaussian 
distribution with mean (Aif^^) = 14 ms and standard deviation 4 ms, and a = 24/ 
(Azi^k). The mean inter-spike interval of such an input approximately equals the mean 
period of the coupling- and input-free neurons (5), and each neuron receives an 
independent random synaptic input of intensity /, which does not significantly 
perturb its natural spiking frequency, see Supplementary Fig. S3. 



Spike timing-dependent plasticity (STDP). The neural ensemble (5) is equipped 
with STDP^^"^^, where the synaptic weights k^ = kyit) are updated in a point process- 
like manner as kij\-^kij + 5-W(^Atij) according to the following STDP function when 
the post-synaptic neuron / or the pre-synaptic neuron j fires: 



W At, 



Aie 



Atii>0 



-A2e'^''i/'\ Atij<0 



(7) 



Atij = ti — tj is the time difference between the nearest spike onsets of the neurons / 
and j. The spike onsets are detected by the upward zero crossing by the membrane 
potential following an inter-spike hyperpolarization. If the time difference Atij is 
negative (the spike of the post-synaptic neuron advances that of the pre-synaptic 
neuron), the corresponding synaptic weight is depressed, otherwise, it is potentiated. 
Since STDP acts on a slower time scale than the underlying neural activity, we 
consider a slow changing rate of the synaptic weights (5 = 0.001, see also Refs. 14, 17. It 
also prevents from spuriously moving away from a branch of attractors known as R- 
tipping'*". The synaptic weights k^ are confined to the interval [k^my ^max] = [0, 0.5] by 
setting kij to /Cmm as soon as ky tends to be depressed below k^^^ via STDP or, 
respectively, to /c^ax if kjj is potentiated beyond this value. 

We consider parameters Ai = 1, A2 = 0.5, Ti = 1.8, and T2 = 6 such that, if the 
input is absent (jJ^p^* = 0), STDP can lead to a multistability of a weakly coupled 
asynchronous state and a strongly coupled synchronized state, see Supplementary 
Section SI A and Refs. 18, 41. Depending on the initial coupling, the mean coupling 
K{t) = J2 ^iji^) saturates at K«0.03 in the weakly coupled regime and at 

K~0.25 in the strongly coupled regime, and the synchronization order parameter^^ 

N 

J2 is on average {R) ~ 0.28 and {R) ~ 0.92, respectively. Here, xjjjit) = 

;=1 

2Ti{t — tj^rn)l{ij,m+\ ~ tj,m)y tj,m ^t< U^m+iy m = 0, 1, 2, . . ., approximates the phase of 
neuron j, and tj^m are the times of its spike onsets. For ij"^^* = 0, the mean coupling 
K'^kmax/2 [Fig. la], and the established uni- directional hierarchical coupling topo- 
logy [Fig. lb] results from the asymmetry of the STDP rule (7) and a phase locking of 
the neurons to each other, see Supplementary Section SIA and Refs. 42, 43. 

Phase oscillators with STDP. We also consider the Kuramoto phase ensemble^ ^ with 
independent noise 



R 



N- 



"^kijsinl^il/j-il/ij+fi^i, 



(8) 



where the coupling weights ky obey the STDP rule (7) with parameters Ai = 1, A2 = 
0.5, Ti = 0.45, T2 = 1.5, and 5 = 0.005, and are confined to the interval ky e [k^my 
kmax\ with /Cmin = 0 and /Cmax = 1- The update (7) of the coupling coefficients is based 
on the phase difference W = W{A\J/y), A\J/y = \J/j — xj/j, between the pre- and 
postsynaptic oscillators; and i, respectively, and undertaken when their phases cross 
zero (mod 2n), which is declared as a spike onset if the phase visited an "inter-spike" 
value n (mod In) just before. The phase-based STDP was chosen in order to avoid any 
influence of the frequency dynamics and to show that the discussed effects are based 
on the noise-induced dynamics of the phase differences. The phase oscillators (8) are 
perturbed by a normally distributed independent noise with {^i) = 0 and {^i{t), 
^jit')} = ^{i ~ ~ t'), and intensity ft. The natural frequencies cOf are uniformly 
and randomly distributed in the interval [0.9, 1 . 1 ] . In calculations for fitting the noise- 
induced distribution of the phase differences we also used a regular distribution co, = 
0.9 + 0.2(/ — l)/(Ar — 1). For numerical integration of model (8) with noise input we 
used the Heun method^"^ with step size 0.01. 

As for the HH neural ensemble, the noise-free (/i = 0) phase ensemble (8) with the 
phase-based STDP (7) demonstrates a multistability of a weakly coupled asyn- 
chronous regime (with K~0.11 and (R) ~ 0.4) and a strongly coupled synchronized 
regime (with K~0.5 and (R) ~ 0.98). In the latter regime a uni-directional hier- 
archical coupling from fast to slow oscillators is established, similarly to the HH 
neural population, see Fig. lb. 

Fixed coupling. To compare the effect of the noise for the HH and phase ensembles 
with STDP to those without STDP, in the latter case the coupling matrix was fixed to 
Kpix = {kij\, where ky = k = k^^Jl for / < j, and ky = 0 otherwise (the oscillators are 
sorted with respect to increasing natural frequency), as for the strongly coupled 
regimes established without noise. 
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